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Abstract 



We give a summary of the recent progress made by the authors and 
. collaborators on the asymptotic analysis of the two matrix model with 

a quartic potential. The paper also contains a list of open problems. 

1 Two matrix model: introduction 

The Hermitian two matrix model is the probability measure 

J_^-nTr{V{M,)+WiM2)~rM,M2) ^^idMs (1.1) 

defined on pairs (Mi,M2) of n x n Hermitian matrices. Here V and W are 
two polynomial potentials, r ^ is a coupling constant, and 

Z„ = / e-nTr(V(Mi)+H/(M2)-rMiM2) ^MidMa 



is a normalization constant in order to make (|1.1|) a probability measure. 

In recent works of the authors and collaborators |20 1 l2 H [22 t l35j the model 
was studied with the aim to gain understanding in the limiting behavior of 
the eigenvalues of Mi as n — )• oo, and to find and describe new types of 
critical behaviors. 
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The results should be compared with the well known results for the 
Hermitian one matrix model 

1 g-nTr(V(M))^^^ (1.2) 

which we briefly summarize here. The eigenvalues of the random matrix M 
from ()1.2p have the explicit joint p.d.f. 

^" i<fc i=i 

which yields that the eigenvalues are a determinantal point process with 
correlation kernel 

n-l 



fc=0 



where {pk,n)k is the sequence of orthonormal polynomials with respect to the 
weight function e~"^(^) on the real line. As n — > oo the empirical eigenvalue 
distributions have an a.s. weak limi10 ^ Z^j=i ~^ where ^* is a non- 
random probability measure that is characterized as the minimizer of the 
energy functional (Coulomb gas picture) 

Evif^) = jj log ^^^^d^x{x)d^l{y) + j V{x)d^x{x) (1.3) 

when taken over all probability measures on the real line. For a polynomial 
V the minimizer /i* is supported on a finite union of intervals [15j. In 
addition there is a polynomial Q of degree deg V — 2 such that 



^z) = V\z) - 
is the solution of a quadratic equation 



diil{s) 



e -V'{z)^ + Q{z) = Q. (1.4) 

From this it follows that /i| has a density with respect to Lebesgue measure 
that is real analytic in the interior of any of the intervals and that can be 
written as 

, , dn*(x) 1 / r-r ^ 

dX TT 



^i.e., for any bounded continuous function /, we have lim — X^'L, fi^j) = f fdfJ.* 
almost surely. 
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where q denotes the negative part of the polynomial 



2 Limiting eigenvalue distribution 

2.1 Vector equilibrium problem 

Guionnet [26j showed that the eigenvalues of the matrices Mi and M2 in the 
two matrix model (jl.ip have a limiting distribution as n — )• 00. The results of 
|26j are in fact valid for a much greater class of random matrix models. The 
limiting distribution is characterized as the minimizer of a certain functional, 
which is however very different from the energy functional (II. 3p for the one 
matrix model. 

Our aim is to develop an analogue of the Coulomb gas picture for the 
eigenvalues of the matrices in the two matrix model (jl.ip . We have been 
successful in doing this for the eigenvalues of Mi in the case of even poly- 
nomial potentials V and W with W of degree 4. Thus our assumptions 
are 

• y is an even polynomial with positive leading coefficient, 



• and r > (without loss of generality). 

We use the following notions from logarithmic potential theory [38] 




which define the mutual logarithmic energy /(^u, v) of two measures and 
V, and the logarithmic energy /(^) of a measure Then the limiting mean 
distribution of the eigenvalues of Mi is characterized by a vector equilibrium 
problem for three measures. This involves an energy functional 



-/(//1,/X2) -/(/X2,//3)+ Vi{x) dfiiix) + / V3{x)dn3ix) (2.1) 



defined on three measures /ii, ^2, f^3- Note that there is an attraction 
between the measures fii and /i2 and between the measures ^2 and fj,^, 



• W{y) 



= iy^ + f y2 with a G M, 






3 



while there is no direct interaction between the measures /ii and ^3. This 
type of interaction is characteristic for a Nikishin system [36]. 

The energy functional (|2.ip depends on the external fields Vi and V3 that 
act on the measures ni and /X3 in (|2.ip . The vector equilibrium problem will 
also have an upper constraint C72 for the measure H2 ■ These input data take 
a very special form that we describe next. 

External field Vi: The external field that acts on /Ui is defined by 

VUx) = V(x) + mm(W(s) - txs), (2.2) 

where we recall that W{s) = ^s^ + ^s^. For the case a = 0, this is simply 
Viix) = V{x) - ||rx|^/3. 

External field V3: The external field that acts on the third measure is 
absent if a > 0, i.e., 

V3(x) = ifa>0. 

The function s G M 1— > iy(s) — txs has a global minimum at s = si{x) 
and this value plays a role in the definition of Vi, see (12. 2p . For q < 0, and 

X £ {—x*{a),x*{a)), where 

2 /_Q,\3/2 

^*(") = 7( — ) ' «<o- 

The function s € M 1-^ — txs has another local minimum at s = S2{x), 

and a local maximum at s = S3(x). 
Then V3 is defined by 

Vsix) = {W{s3ix)) - TXSsix)) - iWis2ix)) - TXS2(x)) (2.3) 
if X G {-x*{a),x*{a)), and 

Vsix) = otherwise. 

Upper constraint CT2: The upper constraint (T2 that acts on the second 
measure is the measure on the imaginary axis with the density 

da2{z) ^ r ^^^^ z£iR. (2.4) 

\dz\ TT s3+as=r2 

In case a = this simplifies to ^ = ^t'^/^\z\^/^ . 
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If a < then the density of (T2 is positive and real analytic on the full 
imaginary axis. If a > then the support of 02 has a gap around 0: 



The following result is Theorems 1.1 in |22] . 

Theorem 1. There is a unique minimizer {^\,^2it^V) of the energy func- 
tional (12. ip subject to the conditions 

(a) ^1 is a measure on M with ^i(M) = 1, 

(b) ^2 is a measure on iM with fi2{i^) = 2/3, 

(c) ^3 is a measure on M with usiM.) = 1/3, 

(d) fJ.2 <cr2, 

with input data Vi, V3, and a as described above. 

The proof of the existence of a minimizer was completed and simplified 
in [27], see subsection 14.21 below. 

Now that we have existence and uniqueness, it is natural to ask about 
further properties of the minimizer. The three measures /i*, o" — /ig and ^3 are 
absolutely continuous with respect to the Lebesgue measure with densities 
that are real analytic in the interior of their supports, except possibly at the 
origin. Furthermore, denoting by S{fi) the support of a measure /i, we have 

• the support of is a finite union of bounded intervals on the real line; 

• there exist C2 > such that S{a2 — ^2) = {—ic2-,ic2)-, and if C2 > 
then the density of 02 — ^1% vanishes like a square root at ±zc2; 

• there exist C3 > such that S{^*^) = R \ (—03,03), and if C3 > then 
the density of /U3 vanishes like a square root at ±03. 

In a generic situation, the density of ^\ is strictly positive in the interior 
of its support and vanishes like a square root at endpoints. In addition strict 
inequality holds in the variational inequality outside the support S{ijl\). 
Moreover, generically if C2 = the density of a — ^2 is positive at the origin, 
and likewise if C3 = the density of /Ug is positive at the origin. If we are 



supp(cr2) = (— ioo, — zy*(a)] U [%*(«), ioo) 



where 




2 



3/2 



a > 0. 
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in such a generic situation, then we say that (V, W, r) is regular. See [22\ 
section 1.5] for more details and a discussion on the singular situations that 
may occur. 

The following is Theorem 1.4 in |22) . 

Theorem 2. Let be the first component of the minimizer in Theorem 1, 
and assume that {V, W, r) is regular, then as n oo with n = (mod3), the 
mean eigenvalue distribution of Mi convergences to . 

We are convinced that the theorem is also valid in the singular cases, 
which correspond to phase transitions in the two matrix model. The condi- 
tion that n is a multiple of three is non-essential as well. It is imposed for 
convenience in the analysis. 

In [22j only the convergence of mean eigenvalue distributions was consid- 
ered, which is a rather weak form of convergence. However, when combined 
with the results of [26] it will actually follow that the empirical eigenvalue 
distributions of Mi tend to fi^ almost surely. 

The analysis of [22j also proves the usual universality results for local 
eigenvalue statistics in Hermitian matrix ensembles, given by the sine kernel 
in the bulk of the spectrum and by the Airy kernel at edge points. In 
non-regular situations one may find Pearcey and Painleve II kernels, while 
in multi-critical cases new kernels may appear. This was indeed proved 
recently in [19], see subsection 14. II below. 

2.2 Riemann surface 

A major ingredient in the asymptotic analysis in |22j is the construction of an 
appropriate Riemann surface (or spectral curve), which plays a role similar 
to the algebraic equation ()1.4p in the one matrix model. The existence of 
such a Riemann surface is implied by the work of Eynard [Mj on the formal 
two matrix model. Our approach is different from the one of Eynard, in that 
we use the vector equilibrium problem to construct the Riemann surface, and 
in a next step we define a meromorphic function on it. 

The main point is that the supports S{fil), S{a — ^'^^ ^(.l^s) as- 
sociated to the minimizer in Theorem 1, determine the cut structure of a 
Riemann surface 

4 

n=\Jn, 
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with four sheets 



7^l = C\5(^D, 
n2 = C\iSifil)US{a2-fi*2)), 
ns = C\{S{a2-fil)US{fil)), 
7^4 = C\5(^^). 

The sheet TZj is glued to the next sheet 7^j+i along the common cut in 
the usual crosswise manner. The meromorphic function on TZ arises in the 
following way, see Proposition 4.8 of |22] . 

Proposition 3. The function 

Uz) = V'{z)- f ^G7^l, 
J z-x 

extends to a meromorphic function on the Riemann surface TZ whose only 
poles are at infinity. There is a pole of order deg V at infinity on the first 
sheet, and a simple pole at the other point at infinity. 

The proof of Proposition [3] follows from the Euler-Lagrange variational 
conditions that are associated with the vector equilibrium problem. See 
section 4.2 of [22j for explicit expressions for the meromorphic continuation 
of ^1 to the other sheets. 

It follows from Proposition [3] that is one of the solutions of a quar- 
tic equation, which is the analogue of the quadratic equation ()1.4p that is 
relevant in the one matrix model. 

3 About the proof 

We describe the main tools that are used in the proof of Theorem [2l 
3.1 Biorthogonal polynomials 

We make use of the integrable structure of the two matrix model that is de- 
scribed in terms of biorthogonal polynomials. In this context the biorthogo- 
nal polynomials are two sequences of monic polynomials {pj,n)j and {qk,n)k 
(depending on n) with degpj^n = j and degqk^n = k, that satisfy 

-oo J — oo 



7 



see [5l [6l O [23l |25] . These polynomials uniquely exist, have real and simple 
zeros [23], and in addition the zeros of pj^n and Pj+i.n interlace, as well as 
those of Qk^n and qk+i,n, see pO] , 

There is an explicit expression for the joint p.d.f. of the eigenvalues of 
Ml and M2 

-^TT^det (^?2 ^ ^i'2/''''^'M (3-1) 

(n!)2 \K^^^'\y,,y^) Kt'\yi,y,)' 



with 4 kernels that are expressed in terms of the biorthogonal polynomials 
and their transformed functions 



00 



J —00 

/oo 

as follows 

K^^''^\xi,X2) = p— Pfc,n(a;i)Qfc,n(a:2), (3.2) 



n-1 



k=0 
n-1 



k,n 



K^n''^\x,y) = ^ T^Pfc,n(a;)gfc,n(y), 



A:=0 

71—1 ^ 

n— 1 ^ 

K^''^\yi,y2) = ^ -To-Pk,n{yi)qk,n{yi)- 

k=0 k,n 

The joint p.d.f. (j3.ip is determinantal, which means that eigenvalue correla- 

(i i) 

tion functions have determinantal expressions with the same kernels Kn , 

i,j = 1,2. In particular, after averaging out the eigenvalues of M2 we get 

that the eigenvalues of Adi are a determinantal point process with kernel 
^(1,1) 

A natural first step to compute the asymptotic behavior of the polyno- 
mials and hence the kernels, is to formulate a Riemann-Hilbert problem (= 
RH problem) for the polynomials. Several different formulations exist in the 
literature [3 ESl EH ES] • The analysis in [2ll ESI [35] is based on the RH 
problem in [32] that we will discuss in the next subsection. 
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3.2 Riemann-Hilbert problem 



It turns out that the kernel (|3.2p has a special structure which relates it to 
multiple orthogonal polynomials and the eigenvalues of Mi (after averaging 
over M2) are an example of a multiple orthogonal polynomial ensemble [3(5]. 
This is due to the following observation of Kuijlaars and McLaughlin [32j . 

Proposition 4. Suppose W is a polynomial of degree r + 1, and let 



k —n(V (x)+W (y)—Txy) 



y e 



dy, k = 0, 



1. 



Then the biorthogonal polynomial pj^n satisfies 



Pj,n{x)J-Wk,n{x) dx = 0, 1 = 0,... 



J 



1, 



(3.3) 



for A; = 0, . . . , r — 1. 

The conditions (|3.3p are known as multiple orthogonality conditions [3], 
and they characterize the biorthogonal polynomials. 

The advantage of the formulation as multiple orthogonality is that these 
polynomials are characterized by a RH problem of size (r + 1) x (r + 1), [40J . 
which we state here for the case r = 3 and for j = n with n a multiple of 
three. Then the RH problem has size 4x4 and it asks for a 4 x 4 matrix 
valued function y on C \ M such that 

(a) y : C \ M ^ C^^^ is analytic, 



1 
1 



(b) Y+ix 



Y^(x) 



/I 




Wo^n(X 
1 




Y^(x) (yi(a;))) denotes the limiting value of Y(z) as z 
upper (lower) half plane, 




VO 



for X G 



where 



X from the 



(c) y(z) = (/4 + o(i)) 
















2— n/3 














^— ra/3 





Vo 











as z 



oo. 



The RH problem has a unique solution which is given by 

fPn,n C{pn,nWo ,n) C (^PnnWl^n 

C{pn}iWo^n) C{pn}iWi^n) C{pn}nW2,n 

'{pnliWi^n) C{pn}iW2,n 

\Pn)n C{pn}iWo,n) C{pn}nWi^n) C{pn}nW2,n)/ 



Y 



JO) 



Pn}n C{pn)nWQ^n) 
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where pn,n is the nth degree biorthogonal polynomial, Pn}i, Pnji, Pn}n are 
three polynomials of degree < n — 1 that satisfy certain multiple orthogonal 
conditions and Cf is the Cauchy transform 

cfiz) =^ r ^dx. 

2m X- z 

By using the the Christoffel-Darboux formula for multiple orthogonal 
polynomials [101 E] the correlation kernel Kn'^^ for the eigenvalues of Mi 
can be expressed in terms of the solution of the RH problem as follows 



K^'^\x,y) = (O ■wo,n(y) ■wi,„(y) W2,n{y)) 



Y-\y)Y+{x) 
2TTi{x — y) 



n\ 




Vo/ 



(3.4) 



Multiple orthogonal polynomials and RH problems are also used for random 
matrices with external source [8l [10] and models of non-intersecting paths 
[31| . In these cases, correlation kernels for the relevant statistical quantities 
are also expressed in terms of the corresponding RH problem through (I3.4p . 

3.3 Steepest descent analysis 

The remaining part of the proof of Theorem[2]is an asymptotic analysis of the 
RH problem via an extension of the Deift-Zhou steepest descent method [16[ 
[T7] . The vector equilibrium problem and the Riemann surface play a crucial 
role in the transformations in this analysis. For the precise transformations 
and the many details that are involved we refer the reader to [22] . Following 
the effect of the transformations on the kernel (13.41). one finds that 



lim -K^^^^\x,x) 



dn\{x) 



>oo n ax 
which is what is needed to establish the theorem. 

A somewhat similar steepest descent analysis is done in [8] for a random 
matrix model with external source, where vector equilibrium problems and 
Riemann surfaces also play an important role. 



4 Further developments 

4.1 Critical behavior in the quadratic/quartic model 

For the case V{x) = ^x"^ the spectral curve can be computed and a classifi- 
cation of all possible cases can be made explicitly. 
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r 



Pearcey 
transition 




-2 -1 



Figure 4.1: Phase diagram for the quadratic case V{x) = ^x^. 

Case I: € S{fil) n S{fil) and S{a2 - /U^); 
Case II: € and S{nl) U 5(^2 - /U^); 

Case III: € 5(f72 - /x*) and S{fil) U 
Case IV: G 5(^5;) and 5(^2 - /i^) U 

Phase transitions between the regular cases represent the critical cases. 

The quadratic/quartic model depends on two parameters, namely the 
coupling constant r and the number a in the quartic potential W{y) = 
^ + a^. Figure UTT] (taken from [20j) shows the phase diagram in the a-T 
plane. Critical behavior takes place on the curves = a + 2 and ar"^ = — 1. 
On the parabola = a + 2 a gap appears around in the support of either 
/i^ (if one moves from Case I to Case II) or /ig (if one moves from Case I to 
Case IV). This is a transition of Painleve II type which also appears in the 
opening of gaps in one matrix models [Oj |i3j . On the curve ar^ = — 1 a gap 
appears in the support of either (if one moves from Case IV to Case III) 
or /ig (if one moves from Case II to Case III), while simultaneously the gap 
in the support of (72 — ^2 closes. This is a transition of Pearcey type, which 
was observed before in the random matrix model with external source and 
in the model of non-intersecting Brownian motions [111 [T2| [39] . 

The phase diagram has a very special point a = —1, t = 1 which is on 
both critical curves, and where all four regular cases come together. For 
these special values, the density of /i* vanishes like a square root at the 
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origin, which is an interior point of >S'(^^). The local analysis at this point 
was done very recently by Duits and Geudens [19]. They found that in 
the asymptotic limit, the local eigenvalue correlation kernels around are 
closely related to the limiting kernels that describe the tacnode behavior for 
non-intersecting Brownian motions [H [T51 [55] . More precisely, the kernels 
can be expressed in terms of an extension of the same 4x4 RH problem in 
[18]. However, they are constructed in a different way out of this 4x4 RH 
problem and, as a result, these kernels are not the same. 



4.2 Vector equilibrium problems 

The analysis in [22] of the vector equilibrium problem was not fully complete, 
since the lower semi-continuity of the energy functional (|2.ip was implicitly 
assumed but not established in [22]. 

In the recent papers [41 127( the vector equilibrium problem was studied in 
a more systematic way, in the more general context of an energy functional 
for n measures 

n n n „ 

E{ni,...,fin) = ^^CijI{ni,Hj) + Y^ / Vj{x)dfij{x), (4.1) 
i=i j=i j=i 

where C = {cij)fj^-^ is a real symmetric positive definite matrix (in [3] 
also semidefinite interaction matrices are considered). The external fields 
Vj : Tij M U {oo} are lower semi-continuous with domains Ej that are 
closed subsets of C. Let mi, . . . , m„ be given positive numbers and assume 
that for every i = 1, . . . ,n, 



liminf | Vi{x) — ( ^ ^CijUij | log(l + |xp) | > — oo. 

Under these assumptions it is shown in |27] that the energy functional (|4.ip , 
restricted to the set of measures with /ij(Ej) = rrij for j = 1, . . . , n, 

(a) has compact sub-level sets {E < a} for every a G M, (so E is in 
particular lower semi-continuous), and 

(b) is strictly convex on the subset where it is finite. 

This guarantees existence and uniqueness of a minimizer of (j4.ip . provided 
that E is not identically infinite. Existence and uniqueness of a minimizer 
readily extends to situations where the domain of E is further restricted by 
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upper constraints /ij < aj for j = 1, . . . ,n, again provided that E is not 
identically infinite on this domain. In particular, this applies to the energy 
functional (j2.ip for the two matrix model with quartic potential with the 
constraint ^2 ^ (^2 described in section 2. 

4.3 Open problems 

There are numerous intriguing questions and open problems that arise out 
of our analysis. 

(a) What is the motivation for the central vector equilibrium problem? 
In the one-matrix model there is a direct way to come from the joint 
eigenvalue probability density to the equilibrium problem. We do not 
have this direct link for the two matrix model. 

(b) How is the vector equilibrium problem related to the variational prob- 
lem from [26J? 

(c) A possibly related question: is there a large deviation principle asso- 
ciated with the vector equilibrium problem? See e.g. |2] for the large 
deviations interpretation of the equilibrium problem for the one matrix 
model. 

(d) Our analysis is restricted to even potentials V and W. This restric- 
tion provides a symmetry of the problem around zero, which is the 
reason why the second measure fi2 in the vector equilibrium problem 
is supported on the imaginary axis. If we remove the symmetry then 
probably we would have to look for a contour that replaces the imag- 
inary axis. It is likely that such a contour would be an S'-curve in 
a certain external field, but at this moment we do not know how to 
handle this situation. See \33\ [37 ] for important recent developments 
around S'-curves for scalar equilibrium problems. 

(e) Extensions to higher degree W is wide open. If deg W = d then one 
would expect a vector equilibrium problem for d — 1 measures. It 
may be that S-curves are needed for d > 6, even in the case of even 
potentials. 

(f) Exploration of further critical phenomenon in the two matrix model. 
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